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Using memory-function formalism we show that in sheared colloidal suspensions the fluctuation- 
dissipation theorem for self-diffusion, i.e Einstein's relation between self-diffusion and mobility ten- 
sors, is violated and propose a new way to measure this violation in Brownian Dynamics simulations. 
We derive mode-coupling expressions for the tagged particle friction tensor and for an effective, 
shear-rate dependent temperature. 
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There has been a lot of interest recently in non- 
equilibrium behavior of colloidal systems pj. On the 
practical side, it has been stimulated by the impor- 
tance of non-equilibrium properties for the preparation 
and processing of colloidal materials. From the more 
fundamental perspective, colloidal suspensions serve as 
model soft-glassy systems: they have properties similar 
to those of more complex soft materials but are simple 
enough to allow for detailed microscopic, experimental 
and theoretical investigations. Additional impetus came 
from analysis of simple statistical mechanical models (i.e. 
fully connected spin systems) that predicted violation of 
the fluctuation-dissipation theorem (FDT) out of equi- 
librium Q, and an intriguing connection between non- 
equilibrium and glassy properties Q . Subsequently, FDT 
violation was found by means of computer simulations in 
model supercooled fluids 0, |]| . 

Although non-equilibrium phenomenology of colloidal 
suspensions is reasonably well described by so-called 
"schematic models" ja, 0, it is of great fundamental 
and practical interest to develop more microscopic ap- 
proaches. First, connection between non-equilibrium and 
glassy behavior, if it also exists for colloidal systems, 
may provide new insights into the glass transition prob- 
lem. Second, simulational and experimental studies of 
colloidal systems provide detailed information that can- 
not be described in terms of schematic models. Third, a 
more microscopic approach would allow one to correlate 
microscopic properties and macroscopic behavior. 

Recent investigations of the non-equilibrium behavior 
can be divided into two categories. Initially, transient be- 
havior of glassy systems, i.e. aging, attracted the most 
attention (j, y|. |_lf. Microscopic, theoretical analysis of 
aging is in its infancy More recentby.relaxation un- 
der steady shear was investigated |a, la IS( ■ Analysis of 
sheared suspensions is easier: stationary nature of the 
shear flow restores time-translational invariance and thus 
simplifies the problem both conceptually and technically. 

Colloidal suspensions under steady shear flow were the 
subject of two recent theoretical investigations [lH lllj . 
Both approaches were based on the "least inadequate" 
[lij l microscopic theory of colloidal dynamics: mode- 



coupling theory (MCT) [13j. The difference between 
them parallels the difference between derivations of MCT 
for non-sheared colloids: Ref. 11(1 used the projection 
operator method whereas Ref. |ll| started from general- 
ized fluctuating hydrodynamics. Both approaches recov- 
ered the most important features of the soft-glassy rhe- 
ology: accelerated relaxation of sheared colloidal fluids 
and shear melting of the colloidal glass. These phenom- 
ena were attributed to flow-induced advection of density 
fluctuations and the resulting perturbation of the "cage 
effect" . Neither work, however, addressed FDT violation. 

The goal of this Letter is to investigate the origin of 
FDT violation for the simplest possible process: self- 
diffusion (i.e. diffusion of a tagged particle) in a sheared 
colloidal suspension [bl . In this case FDT amounts 
to Einstein's relation between the self-diffusion tensor 
and the tagged particle mobility tensor. We general- 
ize the conventional memory function description of self- 
diffusion to colloidal systems under shear and derive a 
Green-Kubo-like relation for the self-diffusion tensor in 
a sheared suspension. Next, we derive a Green-Kubo- 
like relation for the tagged particle mobility tensor. We 
show that FDT violation is associated with the non- 
equilibrium nature of the stationary, shear-rate depen- 
dent probability distribution. We propose a new ap- 
proach to monitor FDT violation in Brownian Dynamics 
simulations that does not require introducing an external 
perturbation. Finally, we use the memory function ap- 
proach to derive MCT expressions for the tagged particle 
friction tensor and an effective temperature. 

We start with the definition of the self-intermediate 
scattering function, F s (ki; k 2 ; t): 



F,(k i; k a ;t) = (n,(ki)exp(nt)n a (-k 2 )> 



(1) 



Here n s (ki) is the Fourier transform of the microscopic 
density of particle number 1, i.e. the tagged particle, 



n s (ki) = e tk i r i. Furthermore, fl is the A^-particle evo- 
lution operator, i.e. the Smoluchowski operator |l6| . 



(2) 
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where Do is the diffusion coefficient of an isolated col- 
loidal particle, (3 = l/(fcsT), F; is a force acting on 
the Zth particle, and v(r) = T ■ r is the shear flow with 
r = 7xy being the velocity gradient tensor (7 denotes 
the shear rate). Finally, in Eq. (• • ■) denotes the sta- 
tionary, shear-rate dependent ensemble average at tem- 
perature T. The probability distribution stands to the 
right of the quantity being averaged and all operators 
act on it as well as on everything else. Note that trans- 
lational symmetry leads to F s (ki; k 2 ; t) oc <Jki(t),k 2 where 
ki(t) = ki +'r • kit with *r being the transpose of T. 

To derive the memory function representation we start 
from the exact expression for the Laplace transform (LT) 
of the time derivative of F s (ki; k 2 ;i), 



LT(F S ) = (n.(ki)n 



n 



n s (-k 2 ) 



(3) 



We define the projection operator on the space spanned 
by the tagged particle density, 



P, 



• n s (-q))(n s (q).. 



(4) 



Note that, in contrast to Ref. [lCj, the definition of P in- 
volves the stationary, shear-rate dependent distribution. 
Next, we follow the usual projection operator manipula- 
tions and arrive at the following memory function repre- 
sentation of the time derivative of F s (ki; k 2 ; t), 

LT(F S ) = ^(n s (k 1 )On s (k 3 ))F s (k 3 ,k 2 ;2) 

k 3 

1 



^(n s (k!)nQ s 



xF s (k 3 ,k 2 ;z), 



Q s fln s (k 3 



(5) 



where Q s is the projection on the subspace orthogonal to 
the tagged particle density, Q s = 1 — P s , 

Using translational invariance of the sheared suspen- 
sion and the fact that in the stationary state the average 
force acting on the tagged particle vanishes, we get 

(n,(ki)fin a (ka)) = -D kf5 kuk3 +k r T~S klM3 , (6) 

oki 

where the second term on the right-hand-side (RHS) of 
Eq. © describes flow-induced advection. 

Furthermore, using the explicit form of the evolution 
operator we obtain the following identity: 



n s (ki)QQ s 



1 



ki-(j.(ki) 



z - Q s ClQi 
1 



z - Q S UQ 



Q s fln s (k 3 ) 

(2jf(-k 3 )-j s (-k 3 )))-k3. 



(7) 



Here j s (k) is a projected tagged particle current density, 
j s (k)=Q J fl (-ik + ^ 1 +v 1 )e- fal , (8) 



and (k) is a projected, effective current density, 

jf (k) = Q S D (-<k + 0Ff) e" ikri . (9) 

In Eq. © Ff ff is the effective force acting on the tagged 
particle that is defined in terms of the stationary, shear- 
rate dependent probability distribution, 



/3Ff = A] ttP *( ri) ..., rjv ) 



(10) 



Combining Eqs. (|6I10|I . in the small wavevector, long 
time limit we get 

LT(F S ) = f -kj D ■ kx +ki • T F s (k l7 k 2 ; z), 

(11) 

where the self-diffusion tensor D is given by the following 
Green-Kubo-like expression: 

/•OC 

D = D Q - ((3D ) 2 / dt (Fx exp(fii) (2Ff - F x )) . 
Jo 

(12) 

Note that in k — > limit projected dynamics (i.e. 
Q s flQ s ) can be replaced by real dynamics (i.e. il) |l7j . 

To get the tagged particle mobility tensor we follow 
approach used by Lekkerkerker and Dhont Il8l . A cal- 
culation along the lines of Sec. Ill of Ref. [18j leads to 
the following Green-Kubo-like formula for the long-time 
tagged particle mobility tensor /j. (here no — (3Dq is the 
mobility of an isolated colloidal particle): 



At = Mo - / dt (Fx exp(fit)Ff ) 



(13) 



Comparison of Eqs. H12fl an d EH shows that Ein- 
stein's relation between self-diffusion and mobility ten- 
sors is violated. The origin of the violation is the differ- 
ence between the force acting on the tagged particle, Fi, 
and the effective force, F\ R . In the absence of the shear 
flow, Fi = F1 S and the usual Einstein relation follows. 
In the presence of the flow, one can follow Ref. [f| and use 
the transverse components of the self-diffusion and mo- 
bility tensors to define an effective, shear-rate dependent 
temperature T off , where fcsT cff — D zz /fi zz . Since we do 
not have an explicit expression for the effective force Ff , 
neither D zz nor /i zz can be obtained from direct simu- 
lational evaluation of its respective Green-Kubo-like ex- 
pression. However, if we obtain D zz from mean-squared 
displacement and measure the force autocorrelation func- 
tion directly, we can obtain the effective temperature: 



k B T cS = 



2D zz k B T 



Do + D zz - ((3D ) 2 J dt (F u exp(fii)F l2 



(14) 

Eq. I|14fl shows that it is possible to monitor FDT vio- 
lation using Brownian Dynamics simulations of the sta- 
tionary, unperturbed state. 
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Before turning to the derivation of MCT expressions 
we first re- write memory function expression Q. We 
define an irreducible evolution operator fl 1 ", 

= -q s (-^ + + v ( r o) 4, (15) 
i 



and then we use standard projection operator manipula- 
tions to obtain the following identity: 



k B T5^ M -£ ^(k 4 ) — (jf (_k 3 )-j s (-k 3 ))^ . (16) 



Note that here (••O -1 denotes the kernel of the inverse 
integral operator; also, £o is the friction coefficient of an 
isolated colloidal particle, £o = l///o- 

Identity (fTT)|) allows us to define the tagged particle 
friction tensor (the inverse of the first factor at the RHS 
of Eq. QIC I) ) and an effective temperature (the second 
factor at the RHS of Eq. (|16(l ). It can be shown that in 
the long-time, small wavevector limit, the former reduces 
to the inverse of the mobility tensor, Eq. (|13fl . 

To derive MCT expressions for the friction tensor and 
the effective temperature we follow the standard proce- 
dure (l^Eof . We project the currents on the part of the 
joint density of the tagged particle and of other particles 
that is orthogonal to the tagged particle density: 

j s (-k) = ^2 "2(-ki,-k 2 )fl'(ki,k2;k 3 ,k4) 

ki,...,k 4 

x(n 2 (k 3 ,k 4 )j s (-k)>, (17) 
jf(-k) « Yl « 2 (-k 1 ,-k 2 ) 5 (k 1 ,k 2 ;k3,k 4 ) 

ki,...,k 4 

x(n 2 (k 3 ,k 4 )jf (-k)>. (18) 

In Eqs. I|17ll8fl ™2(ki,k 2 ) is the part of the joint den- 
sity of the tagged particle and of other particles that is 
orthogonal to the tagged particle density, n 2 (ki,k 2 ) = 
Q s Ez>ie- 4ki ri - lk2 ri , and g(k u k 2 ; k 3) k*) is the in- 
verse of (n 2 (k 3 , k 4 )n 2 (— ks, — ke)). We use fac- 
torization approximation for g, g(ki, k 2 ; k 3 , k 4 ) w 
< 5ki,k3^k2.k 4 (NS(k 2 ))~ , where 5(k 2 ) is the stationary, 



shear-rate dependent structure factor. One should note 
that Eq. (|17|) is exact for pairwise-additive interactions 
whereas Eq. (I18f) constitutes an approximation. 

The average in Eq. I|17f) c & n be expressed in terms of 
the direct correlation force C(k) |21l . 

(n 2 (k!,k 2 )j s (-k)) = -i8 k - klM nD pC(k 2 )S(k 2 ), 

(19) 

(here n is the number density, n = N/V) whereas the 
average in Eq. l|18fl can be expressed in terms the 
non-equilibrium, shear-rate-dependent direct correlation 
function c(k) = (5(k) - l)/(nS(k)), 

<na(ki,k 2 )jf (-k)) = -ik 2 5 k _ kl , k2 n£)oc(k 2 )S , (k 2 ). 

(20) 

Combining Eqs. (|17I20|) with Eq. JTSJ) we can ob- 
tain expressions for the interaction contributions to the 
friction tensor and the effective temperature in terms of 
integrals involving a four-particle, time-dependent cor- 
relation function. We factorize this function in terms of 
the self-intermediate scattering function F s (ki ; k 2 ; t) and 
a collective intermediate scattering function F(ki; k 2 ; t), 

F(k i; k 2 ;i) = i (n(ki) exp(nt)n(-k 2 )) , (21) 

where n(ki) is the Fourier transform of the microscopic 
density, n(ki) = Y^i e ~ As a result we obtain the 

following mode-coupling expressions for the interaction 
contributions to the friction tensor and the effective tem- 
perature in the long-time (i.e. z — > 0) limit: 



•/?£g(j fl (ki)(n irr ) -1 jf (-k 2 )) « - Yl J *4 1 -k3,k 4 C(k4)F s (k 3 ;k 5 ;t)F(k4;k 6 ;t)c(k6)k 6 5 k2 _ k5lk6 , (22) 

k 3 ,...,k 6 

eo(j,(k 1 )(r! iri y 1 (jf(-k 2 )-j 1 (-k 2 ))\ « V / O °da kl -k3,k 4 C(k 4 )^(k 3 ;k 5 ;i)F(k 4 ;k 6 ;0 

N ' v ,.Jo 



k 3 ,...,k 6 
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x (c(k 6 )k 6 - /?C(k 6 )) <5 k2 -k 5 ,k 6 
I 



(23) 



Expression (1221 differs from one derived before jl4| : 
one of the vertices in Eq. i]22li involves the direct correla- 
tion force, whereas the other involves the non-equilibrium 
direct correlation function. In the expression obtained in 
Ref. |LJ] both vertices were identical and given by the 
equilibrium direct correlation function. One should note 
that in equilibrium C e9 (k) = ksT]s.c eq (k); hence the pre- 
vious work implicitly used equilibrium approximation for 
the vertices. 

Expression 12/it shows that FDT violation is associ- 
ated with the non-equilibrium character of the stationary, 
sheared state. In particular, if equilibrium approximation 
for the vertices is used, no FDT violation is obtained. 

One of the most interesting simulational findings is 
that below the MCT transition temperature FDT vio- 
lation persists in the limit of the vanishing shear rate 
0. This can be qualitatively understood on the basis 
of Eq. i|23|) : in the 7 ^ limit the second vertex in 
(12.'it vanishes; however, in the same limit characteristic 
relaxation times of the self and collective intermediate 
scattering functions diverge. Thus, it is possible that 
the expression l|23[) reaches a finite limit as 7 — > 0. In 
order to prove that this indeed happens one needs to cal- 
culate the second vertex in Eq. 1)23(1 . To this end, it 
might be possible to use an approach proposed by Fuchs 
and Cates [ljj: using mode-coupling theory to calculate 
steady state properties by starting from the equilibrium 
state and considering transient dynamics. 

Another very interesting result of Ref. [j| is that in the 
7^0 limit the same effective temperature is obtained for 
different wavevectors and even for different observables. 
Expression l|23|) . in principle, allows one to verify this 
fact theoretically. In particular one could check whether 
tensorial quantity Q23J) reduces to a scalar one in the 
7 — > limit. Again, in order to investigate this, the 
second vertex is needed. 

To summarize, we have derived Green-Kubo-like for- 
mulae for the self-diffusion and mobility coefficients, and 
mode-coupling expressions for the friction tensor and the 
effective temperature. The numerical analysis of these 
expressions is left for future work. 
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